Proportion surviving after 5 years | stcox/stcurve

[Girish Venkataraman]

This seems like a simple question. But I am confused how to do it. I did a 3 parameter Cox model adjusting for two covariates using stcox and then did stcurve at 2 values of my covariate of interest (rt1chop of 0 or 1).

I want to report the proportion surviving at end of 5 years for the two group. What is the best way to get the percentage? I did one below. Can share the stsetting variables if that is needed

stcox ib0.rt1chop i.ecog2 age65, noshow


stcurve, surv at1(rt1chop = 0) at2(rt1chop = 1) ///
lcolor(navy dkorange) lpat(solid -) lwidth(thin thin) ///
graphregion(margin(vsmall) icolor(white) color(white) ) ///
title("CHOP-like vs. EPOCH-like for RT first line") bcolor(black) ///
subtitle("After adjusted P-H regression") ///
note("Adjusted for age and ECOG PS") ///
xtitle("Time since transformation (years)", size(small)) ///
xlabel(0(1)5, labs(small)) ylabel(0(0.25)1, angle(0) labs(small)) ///
ytitle("Overall survival", size(small) margin(medium)) ///
legend(order(1 "CHOP-like" 2 "Intensive regimens") ///
col(1) size(small) symxsize(*0.5) margin(small) rowg(*0.1) ring(0) pos(1)) ///
name(stcurve_byrt1chop, replace)

* Example generated by -dataex-. For more info, type help dataex
clear
input byte(rt1chop ecog2 age65 _st _d) int _origin byte _t0
. . 0 1 1 20619 0
0 0 1 1 1 20626 0
1 0 0 1 1 19876 0
. 0 1 1 0 22302 0
. 0 1 1 1 22186 0
. 0 1 1 0 22111 0
0 1 1 1 1 22335 0
. 1 1 1 1 22179 0
. . 0 1 1 22137 0
1 0 1 1 1 21880 0
. 1 1 1 1 21011 0
. 0 1 1 1 22286 0
0 0 0 1 1 21816 0
0 0 1 1 1 22175 0
1 0 0 1 1 21614 0
0 0 0 1 0 20370 0
. 0 0 1 0 21976 0
0 . 1 1 1 20879 0
1 0 0 1 1 21836 0
0 1 0 1 1 21434 0
1 0 0 1 0 22518 0
0 0 1 1 0 21700 0
. 0 0 1 1 20230 0
. 0 0 1 0 19302 0
1 0 1 1 1 21544 0
0 0 1 1 0 21817 0
. . 0 1 0 20914 0
1 0 1 1 1 21635 0
. 0 0 1 1 20783 0
0 0 0 1 0 22148 0
. 0 0 1 0 22200 0
. 1 1 1 1 20977 0
0 0 1 1 0 22571 0
1 0 0 1 1 19269 0
. 0 0 1 0 21424 0
0 . 0 1 1 19000 0
. 0 1 1 0 22133 0
. 1 1 1 0 21258 0
. . 0 1 0 20800 0
1 . 1 1 0 22432 0
0 0 1 1 1 21069 0
0 1 1 1 1 21084 0
1 1 0 1 1 20237 0
. 0 1 1 1 21627 0
0 0 0 1 1 19757 0
. 0 0 1 0 22137 0
. 0 0 1 0 22553 0
0 0 1 1 0 21097 0
0 1 1 1 1 21833 0
. . 0 1 0 21140 0
0 0 0 1 1 21971 0
. 0 0 1 0 19368 0
. 0 1 1 0 21864 0
0 1 1 1 0 20620 0
0 0 0 1 0 22595 0
. 0 1 1 1 22391 0
1 0 1 1 1 19170 0
1 0 1 1 0 20173 0
. . 0 1 0 22040 0
0 0 1 1 1 21895 0
1 0 1 1 1 22354 0
. . 0 1 0 22715 0
. 0 1 1 1 20251 0
. 0 1 1 0 22102 0
. 0 0 1 0 19606 0
. 0 1 1 1 19067 0
1 . 1 1 0 20236 0
. 0 1 1 0 20200 0
. 0 0 1 0 19586 0
. 0 1 1 1 21536 0
0 0 1 1 1 20804 0
. . 0 1 1 17106 0
. 0 0 1 1 17269 0
0 0 1 1 1 20983 0
0 0 1 1 0 20660 0
. . 0 1 1 16919 0
. . 0 1 0 15937 0
0 1 1 1 1 22082 0
. . 0 1 0 15013 0
0 0 1 1 1 20838 0
0 . 0 1 1 21159 0
0 0 0 1 0 20248 0
. . 1 1 0 17722 0
. 0 1 1 1 21726 0
. 0 1 1 1 20009 0
0 . 1 1 1 18107 0
. 0 1 1 0 20895 0
0 0 0 1 1 19429 0
0 . 0 1 1 19087 0
0 0 1 1 0 19359 0
1 0 0 1 1 20306 0
0 0 0 1 0 17729 0
. 0 0 1 0 17889 0
0 0 1 1 1 21224 0
. 0 1 1 0 20886 0
. . 1 1 1 21964 0
. . 0 1 1 18382 0
. . 1 1 1 17176 0
. . 0 1 0 15767 0
. 0 1 1 1 21063 0
end
label values age65 age2
label def age2 0 "<65 years", modify
label def age2 1 ">=65 years", modify